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We derive the leading one-loop contribution to the one-pion exchange and short-range two-nucleon 
^^ ■ electromagnetic current operator in the framework of chiral effective field theory. The derivation is 

ITi I carried out using the method of unitary transformation. Explicit results for the current and charge 

■ densities are given in momentum and coordinate space. 

|JZ PACS numbers: 13.75.Cs,21.30.-x 

' There have been several recent studies on the nuclear exchange electromagnetic currents within the framework of 
chiral effective field theory pHil, see also Q for an older calculations which, however, is limited to the near-threshold 
kinematics. These studies constitute a natural extension to photon-induced reactions of the theoretical framework 
*^ , formulated by Weinberg two decades ago [y], see [3] for a recent review article. To derive the exchange currents 
^ii^ ■ from the most general effective chiral Lagrangian the authors of Refs. [ll 0, |j| used the framework of "old-fashioned" 
vQ I time-ordered perturbation theory along the lines of [6|. This approach leads, in general, to explicitly energy-dependent 
f~^ ■ potentials and currents. Such energy dependence might cause difficulties in few-body applications. To obtain energy- 
independent nuclear potentials we employed in Refs. [8|, |9| the method of unitary transformation. In Ref. [3|, we 
applied this approach to the long-range parts of the leading two-pion exchange contributions to the current and 
charge densities. In this manuscript, we derive all remaining contributions to the two-nucleon current and charge 
densities at the leading loop order (i.e. of order eQ with Q ~ M^^ referring to low external momenta). 
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I. INTRODUCTION 



.5^ ] It is important to emphasize conceptual differences between our work and the one by Pastore et al. [llS, S|- These 

S^ ■ authors limit themselves to deriving the momentum dependence of the one-pion exchange current and charge operators 

\^ ' at the leading loop level in chiral effective field theory without considering renormalization. Consequently, the values 

. 5t 1 of the various low-energy constants (LECs) entering their expressions cannot be taken from other sources such as 

e.g. pion-nucleon scattering. One, therefore, looses one of the greatest strengths of the effective field theory approach, 

namely the ability to relate different processes. The calculation presented in our work is more ambitious aiming at 

the derivation of renormalized expressions for the exchange current and charge operators. This is a highly nontrivial 

task for the one-pion exchange contributions. Contrary to the calculations in the Goldstone boson and single-nucleon 

sectors, one is dealing here only with an irreducible part of the amplitude (giving rise to nuclear forces and currents) 

which itself is not an observable quantity and is affected by unitary transformations. On the other hand, there is no 

freedom in absorbing the divergences generated by the loop corrections to the one-pion exchange operators since all 

/3-functions of the corresponding LECs are fixed and well known. As we will demonstrate in this work, it is indeed 

possible to exploit the above mentioned unitary ambiguity in such a way that all divergences emerging from pion 
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loops are indeed absorbed by redefinition of the LECs leading to the finite result for the current and charge operators, 
where the values of renormalized LECs can be taken from other sources. 

Our manuscript is organized as follows. In section |TT] we provide a short summary of the method of unitary transfor- 
mation (UT) and explain very briefly the adopted power counting scheme. The effective Lagrangian employed in our 
calculation is specified in section lllll The results for various contributions to the one-pion exchange current and charge 
densities are discussed in detail in section HVl Section |V] deals with the derivation of the short-range contributions. 
A comparison between our work and the calculations by Pastore et al. is presented in section I VII The results of 
our work are summarized in section fVIII The expressions for the relevant terms in the effective pion-nucleon-photon 
Hamiltonian density are listed in appendix|Al while appendix [B] collects the expressions for the relevant loop integrals. 
The Expression for the current and charge density in configuration space are given in appendix [Cj 



II. ANATOMY OF THE CALCULATION 

The derivation of the electromagnetic nuclear current operators is carried out along the lines of Ref . Q , see also [l3| . 
The main steps are summarized below. 

• We begin with the effective chiral Lagrangian in the heavy-baryon formulation and express it in terms of 
renormalized pion and nucleon fields and apply the canonical formalism along the lines of Ref. |11| to derive the 
corresponding Hamilton density. The contributions from tadpole diagrams are taken into account by performing 
normal ordering of the resulting Hamilton density. Notice that the terms in the effective Lagrangian/Hamiltonian 
involving two and more insertions of an external electromagnetic field A'^ are not taken into account since we 
restrict ourselves to the one-photon-approximation (however, the method can straightforwardly be generalized to 
two-photon processes such as Compton scattering off light nuclei) . The obtained contributions to the Hamilton 
density are listed in appendix lAl 



• 



To decouple the purely nucleonic subspace of the Fock space from the rest we apply an appropriately chosen 
UT 



^--'™-(1Ma) 



(2.1) 



Here, 77 (A) denote projection operators onto the purely nucleonic (the remaining) part of the Fock space 
satisfying 77^ = 77, A^ = A, 77A = Xij = and A -I- 77 = 1. The resulting nuclear Hamiltonian rjHrj gives rise to 
the chiral potentials in Refs. [3, [ij, [23, [23. Both the UT and the transformed Hamiltonian are calculated by 
making a perturbative expansion in powers of Q/A, with Q and A referring to the soft and hard scales of the 
order of the pion and p-meson masses, respectively. The power counting is most easily formulated in terms of 
the canonical field dimension k of the interaction vertices. 



00 „ 



• 



K=l 

Here, di, rii and Pi refer to the number of derivatives or M^-insertions, nucleon field and pion field operators, 
respectively. The explicit form of the strong part of the unitary operator [/, i.e. the one in the absence of the 
external electromagnetic field, sufficient to derive the nuclear force up to N^LO is given in Ref. [13 |. 

The effective nuclear current operator r]J'^{x)ri acting in the purely nucleonic subspace of the Fock space is 
defined according to [3, [2^ 

J^{x)^T^U^J^^Jx)Ur,. (2.3) 

Here, Jj^^^^ix) denotes the hadronic current density which enters the effective Lagrangian C^^N-y describing the 
interaction of pious and nucleons with an external electromagnetic field A^. It is given by 

^ba..c(^)-«'^5(5^^^) Q^^ ■ (^-4) 



The A-components of the effective current operator do not need to be taken into account as long as one stays 
below the pion production threshold. The above definition of r]J^{x)r] does, in fact, not fully incorporate the 
freedom in the choice of UT. In particular, one can introduce 77-space UTs 7/[/'ry that depend explicitly on the 
external electromagnetic field Afj, such that 



rjU'r] 
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A^=0 



(2.5) 



Applying such UTs on the nuclear Hamiltonian rjHrj will generate further contributions to the nuclear current 
operator. The resulting ambiguity is analogous to the one in the strong sector which is described in detail in 
Ref. [I2I |l3| . As will be shown below, renormalizability of the one-pion exchange contributions at the one- loop 
level strongly restricts the ambiguity in the definition of 7]J^^{x)rj. 

• The final step in the derivation involves evaluating the emerging loop integrals and expressing the current 
operator in terms of renormalized low-energy constants. This is carried out within the framework of dimensional 
regularization (DR) which allows us to adopt the known expressions for the /3-functions of the LECs entering 
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In the following sections, the various steps in the derivation of the current will be discussed in detail. 



III. EFFECTIVE LAGRANGIAN 



In this work we employ the standard heavy-baryon formulation for the effective Lagrangian. The terms needed in the 
calculation of the leading loop corrections to the one-pion exchange and short-range current operator read |14l - [2(]| 
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+ e^^p^v^fP [Li {N^S^T^N^N^N, - N.S^N^N.t^N,) + L2N,S" N,N,N,] + . . . , (3.1) 

where v denotes the nucleon four-velocity, ( ) stands for the trace in the flavor space and the spin vector is defined as 

% % 1 

Sf, ^ -^'ybCTf^-yv" , cTp^ ^ -[jf„j^], {Sf,,S^} ^ -{Vf^v,,- gf,^) , [Sf„S^]^iep^p„vPS'' , (3.2) 

with the last two relations holding in four dimensions. Further, F, rh and cja refer to the pion decay constant, nucleon 
mass and the nucleon axial-vector coupling in the chiral limit while h, di, Cs,t, oh and Li^2 are further LECs. Notice 
that we only list those terms in the effective Lagrangian which are explicitly needed in our calculations. For example, 
we omit all terms in C\^jq proportional to the LECs Ci as they lead to vertices with at least two pions^ and thus will 

not contribute to the current operator up to the leading-loop order. We further emphasize that the terms in C)^^ 
do not correspond to the minimal set, see [20, [2l| for more details and relations between the different Ui. We will 
address this issue and list the minimal set of contact interactions the nucleon rest-frame at the end of this section. The 
superscript i in Cj^\q, C-nl and Ly^q refers to the number of derivatives and/or quark mass insertions. The unitary 
2x2 matrix U parametrizes the Goldstone Boson fields and is given by 

where ^ is a constant representing the freedom in the definition of the pion fields. The popular cr-model gauge and 
exponential parametrization of the matrix U correspond to ^ = and ^ = 1/6, respectively. Notice that physical 
observables calculated using the effective Lagrangian are, clearly, independent on a particular parametrization of U . 
The quantity x+ is defined via 

X+ = u^xu^ + ux^u , X = '^BM = A/^l2 (3.4) 

with B and M = diag(TO„, md) being a constant and the light quark matrix accounts for the explicit chiral symmetry 
breaking and gives rise to the pion mass 

M^ = M'^{1 + 0{M^)) . (3.5) 

The covariant derivatives of the pion and nucleon fields are defined by 



u^, = i[u^ (dp -irp)u-u{dp~i£p)u^ , 



DpN, = 
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Tp = - [u\df, - irf,)u + u{dp - iif,)u''] , (3.6) 



where r^, Z^ and v^ denote the external right-, left-handed and isoscalar vector currents, respectively, and u = yU. 
The derivative operators Dp and Dp entering £j^^ are defined via 

N^l^pN, = N, (dpN,) + N, (Vp - w(,^)) N, , nJ^pN, = {dpN,) iY„ - N, {Vp - iv^'^^ N, . (3.7) 

Further, f]^^ and Vpy denote the field strength tensors associated with external left-, right-handed and the isoscalar 
currents. 

If. = d^r. - duTp - I {Tp, r,\ , /^^^ = dpi, ~ djp - i [Ip, I,] , v^;^ = dpvl^^ - d^v^:^ , (3.8) 



^ The only exception is the ci-term which also has a contribution that does not involve pion field operators. This contribution can be 
absorbed into redefinition of the nucleon mass. 



while the corresponding covariantly transforming quantities /r^ which enter the pion-nucleon Lagrangian are defined 
according to 



ft = "^ (/"- + ^i- ) " ± " (/m- + «Mt^ 
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We also used traceless matrices /±, defined according to /±^ = f±;^, — {f±^)/2. In this work, we are interested in 
describing the coupling to an external electromagnetic field. In that case, the left- and right-handed currents r^ and 

(s) 

l^ and the isoscalar current vjj, have to be chosen as 
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where Af_i refers to the electromagnetic four-potential. 

We now turn our attention to the Lagrangians Cj^'j^f involving four nucleon field operators. At the order considered, 
there is no need to account for terms involving pion fields. Notice further that Poincare covariance implies that only 7 
out of 18 constants ai are independent and, in addition, also determines the coefficients in front of the leading 1/m^ 
corrections to contact terms, see 0, [20, l2l[ for more details and explicit expressions. In the power counting scheme 
we adopt in the present work, the nucleon mass is treated as a heavier scale compared to the breakdown scale of 
the chiral expansion, see Refs. [3, Q for more details. Accordingly, there is no need to take into account the leading 
relativistic l/m^-corrections to the short-range two-nucleon current at the order we are working. Switching to the 
rest-frame of the nucleon with w^ — (1,0,0,0), making use of the partial integrations and incorporating constraints 
due to the Galilean invariance allows to express the Lagrangian for contact interactions in the standard basis in terms 
of Ci,...,7 used e.g. in [1,1, [l^: 
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where we have introduced 

e=i-t^, A^^Va^ A^^VA-A^t^A^. 
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Notice that we only kept terms at most linear in the electromagnetic four-potential. 
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As already pointed out in the previous section, the derivation of the exchange current operator is carried out using the 
method of unitary transformation which requires the knowledge of the Hamilton density and the Noether currents. 
The transition from the Lagrangian to the Hamiltonian is achieved employing the standard canonical formalism. An 
extended discussion on this can be found in Refs. [3,li3, llll- AH terms in the resulting Hamilton density which enter 
the calculation are listed in appendix 1X1 



IV. ONE-PION EXCHANGE CURRENT 



Wc now turn to the derivation of the two-nucleon electromagnetic current due to a single pion exchange. In section 
IIV Al the derivation and explicit results for the leading loop contributions are presented. Tree-level contributions 
and the renormalization are considered in sections llV Bl and IIV CI respectively. Next, in section HVDI we discuss the 
leading relativistic corrections. Final results for the one-pion exchange current and charge density both in momentum 
and coordinate spaces are summarized in section HV El 



A. Loop contributions 



Following Ref. Q , we classify various loop contributions according to the powers of the LEG gA and the type of the 
hadronic current J;^q, J^j^ or Jq2 as shown in Fig.[TJ Here and in what follows, we adopt the notation of Refs. |j,[l2|. In 

particular, the subscripts a and h in H^ and Jj^f, refer to the number of the nucleon and pion fields, respectively, 
while the superscript k gives the dimension of the operator defined in Eq. (|2.2p . Notice that while class-3 terms 
proportional to g'^ and involving an insertion of Jq2 contribute to the two-pion exchange current, they do not generate 
one-pion exchange diagrams. On the other hand we now have additional contributions from class-8j9 terms which 
do not contribute to two-pion exchange diagrams and, for that reason, were not considered in Ref. [3[. We further 
emphasize that, strictly, speaking (i.e. according to the power oi gA), these diagrams belong to class 5. The algebraic 
structure of the current operator in Fock space in terms of Hat and J^j, for seven classes is given in appendix A of 
Ref. Q. The new terms corresponding to the classes 8 and 9 have the form: 
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Here, the superscript i of A' refers to the number of pions in the corresponding intermediate state. Further, E^ denotes 

the total energy of pions in the corresponding state, E^ = Xii V ^* ^ ~'~ "^^' ^i^h li the corresponding pion momenta. 
We remind the reader that the representation for the power counting in terms of the canonical dimension k allows 
one to easily read off the chiral order associated to a given contribution by simply adding together the dimensions k 
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As already pointed out in section |TT] and in Ref. [3| , we have to employ additional UTs in the ry-space in order to 
maintain renormalizability of the one-pion exchange contributions, see Refs. 12] for a related discussion. For the case 
at hand, one can distinguish between the strong UTs and the ones depending on the electromagnetic four-potential 
A. The general form of the strong UTs up to the considered order in the chiral expansion is given in Ref. [l3]- These 
continuous UTs are parametrized in terms of some (a-priori arbitrary) "angles" di.'^ These parameters turn out to 
be strongly constrained if one requires that matrix elements of the resulting nuclear potentials can be made finite by 
means of redefinition of certain LECs, i.e. if one demands renormalizability at the level of the nuclear Hamiltonian. 
For the UTs considered in Ref. [Q, this condition was shown to lead to a unique expression for the four- nucleon 
force which does not depend on di any more. Similarly, the expressions for the two-pion exchange current operator 
obtained in Ref. 3] are also cii-independent. The additional electromagnetic UTs have not been discussed in that 
reference as they turned out not to affect the two-pion exchange contributions. As will be shown below, it is necessary 
to employ such additional UTs to maintain renormalizability of the one-pion exchange current. To be specific, we 
consider the 77-space UT of the form 
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In that reference, the angles were denoted by a^. 
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FIG. 1: Leading loop contributions to the one-pion exchange current operator. Sohd and dashed lines refer to nucleons and 
pions, respectively. Solid dots are the lowest-order vertices from the effective Lagrangian while the circle-crosses represent 
insertions of the electromagnetic vertices as explained in the text. Diagrams resulting from interchanging the nucleon lines are 
not shown. 

where S is an anti-hermitian operator acting in the 77-space, S = rjSri, S^ = —S. At the order considered, this operator 
can be parametrized as 



5 = ^ft^, 
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with /3i being arbitrary constants and 
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The action of these UTs onto the one-pion exchange contribution to the lowest-order effective Haniihon operator, 
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with -^20 denoting the nonrelativistic kinetic energy term, induces additional, /3-dependent class-2, class-5, class-6 
and class- 7 contributions: 
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ft turns out to be convenient to express /34,5,6 in terms of another set of constants /?, 7 and 6 defined as: 



(4.6) 



(4.7) 



Already at this stage we emphasize that three of the seven parameters, namely 7, /32 and /?7, do not affect the leading 
one-loop contributions to the one-pion exchange matrix elements. 

After these preliminary remarks, we are now in the position to discuss the results for the one-loop contributions. Here 
and in what follows, the expressions for a class-X contribution J^ refer to the matrix element defined according to 



{P1'P2'\J''\P1P2) - 5{pi'+P2'-pi-p2 - k) [Je\ + (1 ^ 2)] 



(4.8) 



Here and in what follows, Pi {pi') refers to the initial (final) momentum of the nucleon i. We will also frequently 
use the momentum transfer variables qi,2 = Pi, 2' — Pi,2- The expressions for the two-pion exchange current and 



charge densities were given in Ref. Q in terms of the most general set of spin-momentum vector and scalar operators 
O1...24 and of .8 as well as isospin operators T1...5. We found that this representation leads to unnecessarily involved 
expressions in the case of the one-pion exchange and short-range currents. We, therefore, refrain from using the 
operators O1...24, Of g and T1..5 in the present work. 

Evaluating matrix elements of the operators in the Fock space as discussed above, we obtain the following results for 
the matrix elements of the current density: 
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'32F4 
and the charge density: 



d^l - l-{qi- 92) 



(7? -f A/2 g| + Af 2 7 (2^)3 ^^0J.iuj++LU-) 



(4.9) 
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C ffA ^3 ^2 • 92 



8F4 ^ g2 _,. ^J2 y (2^)3 V 
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(4.10) 



where 



-I 



{l±k\ 



-AMi 



LO, 



P + Ml . 



(4.11) 



The class-3,4,6 contributions to the current density and class-1, 2,3,4, 5, 8,9 contributions to the charge density are 
found to vanish. 



B. Tree level contributions 



The loop contributions considered in the previous section do not involve the ones emerging from pion tadpole dia- 
grams. These must be explicitly taken into account if one wants to use the values of the renormalized LECs such 
as di determined from e.g. the pion-nucleon system. The treatment of the pion tadpoles in the method of unitary 
transformation is discussed in detail in Ref. 10]. The pion tadpole contributions emerge from contractions of the pion 
field operators when performing the normal ordering of the effective pion-nucleon Hamiltonian and simply lead to 



10 



+ 



+ 



+ 



FIG. 2: Schematic representation of the renormalization of the pion field and the operators -^21 i "^21 '^^'^ "^02 

additional vertex corrections. Following Ref. [10| . we work with renormalized pion field and mass defined according 
to 



TT^ = Z^ ' TTa , Ztt = I + SZ-;^ , 



Ml 



M^ + SMl , 



(4.12) 



where a denotes the isospin quantum number and SZ-^, SM^/M^ ~ 0{Q^ /A^). At the leading loop order, SZ-^ and 
5M^ are given by [l3| 



5Z^ 



Ml = M^ 



F2 
1 



hMl , l-8g 
2F^ 



F2 



A^ 



where the quantity A^ is defined in Eq. (|B.5[) . Notice that in Ref. [Ifll we used the parametrization of the matrix 
U with ^ = 0. We further emphasize that there are no pion self-energy diagrams since we work with renormalized 
pion fields. All effects due to pion self-energy and/or tadpoles are taken into account by vertex corrections in the 
normal-ordered effective Hamiltonian. This is schematically visualized in Fig. [5] More precisely, replacing tTq — > tt^ 

and M^ — > M"^ in £i,f and £^^ generates corrections to Cir-l and C\}q (and, of course, in the corresponding Hamilton 
densities) driven by SZ-^ and SM^. Further corrections, JnOj to the operators H21 , J21 and Jq2 emerge from taking 
normal ordering on the operators -^23 i "^23 ^^"^ "^04 Together with the the wave function renormalization of the pion, 
we obtain the following shifts: 
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2p2^- 



(4.13) 



We point out that, as expected, none of the renormalized operators depends on the (arbitrary) value of ^. 



After these preliminary remarks, we are now in the position to discuss the tree-level contributions to the one-pion 
exchange current and charge densities. The formal operator structure is given by 
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(4.14) 
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FIG. 3: Contributions of the counter terms. Solid dots are the lowest-order vertices from the effective Lagrangian while the 
crosses represent insertions of the electromagnetic vertices as explained in the text. For remaining notation see Fig. [T] 

These operators give rise to diagrams shown in Fig. |31 The explicit form of all vertices entering this expression can 
be found in appendix |^ Evaluating the corresponding matrix elements we obtain the following expressions for the 
current density 



Jtree = ^6— :^ (^81"! + dg (fi • T2)) 



m X q2\ 



^22 fc X [qi X O^l] + <7l 

ai ■ qi 



qi 
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2d2i k X [q2 X cti] 



(4.15) 



while the contributions to the charge density vanish. This is consistent with the fact that the loop contributions to 
the charge density do not contain logarithmic ultraviolet divergences. 



Renormalization 



The expressions given in the previous sections are written in terms of bare parameters and contain ultraviolet- 
divergent pieces. These divergences are cancelled after expressing the bare parameters M, c/a, F, U and di in terms 
of the corresponding renormalized quantities. When carrying out renormalization, one should also take into account 
the contribution induced by the leading-order {O (eQ"^)) one-pion exchange current shown in Fig.|4] 



J 



(eQ ') Wa r- ^ -13 '?2 -92 



^ en • 91 
<ii- 



0-1 



(4.16) 



qi+Ml 
when expressing the ratio (jaI F in terms of the physical LECs qaJF-^. The chiral expansion of this ratio has the form 

Ml, Ml 



9 A _ 9 A 

F^~ F 



F2 



£f = 2^ ( 1 - r^A, - '-^h + ^'-^du 



(4.17) 



F^ 9A 

Clearly, this relation holds modulo higher-order corrections. The resulting induced correction at order O (eQ) reads: 

Ml , . Ml 



ji^Q) ^ e^ ^^^ [fi X r2f 
11 ' 



9172 
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Ml 
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F2 



9A 



(4.18) 



Notice that at the order considered, one can safely replace (ja and F by the corresponding renormalized quantities in 
all expressions given in sections IIV Al and IIVBI 
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FIG. 4: Lowest-order contributions to the pion exchange current operator: the pion-in-flight and seagull graphs. For notation 
see Fig. [T] 

Consider now the LECs li and di which can be decomposed into the divergent parts and finite pieces as follows: 



1 



1 



h = '[(m) + ItL = Y^h + ItL + 7*Y^ log ( — ) ' 






A 
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d,. = <(m) + ^L - d, + j^L + ^^^2p2 log 1 — 






(4.19) 



where the divergent quantity L is defined in Eq. (IB.6| . The corresponding coefficients Pi and 7^ in the framework of 
dimensional regularization (DR) are well known [ij, ll^j Il3) 13 and read: 



^8 = /39 = /3i8 = /322 = 0, /3i6 = -5A + g% hi = -5L H = 2, 



76 



(4.20) 



The expressions for all loop integrals that enter the calculation in DR can be found in Appendix [B] The only exception 
is the part of the class-7 current proportional to the constant 5, for which we did not succeed to find a closed expression. 
Inserting the DR expressions for the integrals entering Eqs. (14.91) . (I4.10|) and the pion tadpole contributions discussed 
above and replacing the bare LECs in terms of renormalized ones, one observes that indeed almost all divergences 
cancel. The only remaining divergent part of the current reads 
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(4.21) 



This implies that we have to choose (3i — \ and — 2 + 2/3 + (5 = in order to be able to renormalize the current 
operator. Here and in what follows, we adopt the choice (5 = and /? = 1. 



D. Relativistic corrections 



Last but not least, we now discuss the leading relativistic corrections. These emerge from the operators in Eq. (I4.14p 

with the vertices ii/21 > ^^20 and J21 being replaced by the corresponding relativistic corrections i^2i ' "^20 ^'i'^ "^21 i 
respectively, whose explicit form is given in appendix [Xj In addition, there are contributions emerging from insertions 

~ (2) 
of the kinetic energy of the nucleon -^20 which have the form 
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(4.22) 



where the constants /3i are defined in Eq. (I4.3p . The additional ry-space UTs considered so far did not involve l/niN- 
corrections. The unitary ambiguity of the leading relativistic corrections can be parametrized in terms of the following 
two additional UTs: 



S = PsSs + PgSg , 

with two new constants /Sg and Pg and the operators 6*8.9 given by 



(4.23) 
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(4.24) 



Notice that the operator Sg with -^21 being replaced by ^21 vanishes which is why the corresponding UT was not 
considered in section IIVBI The effects of these UTs in connection with the nuclear potentials and currents have 
already been investigated, see [2J, |2J| and references therein. In particular, these UTs afi'ect l/?Ti^-corrections to 
the one-pion exchange and 1/rMAr-corrections to the two-pion exchange nucleon-nucleon potentials which appear at 
N'^LO in the chiral expansion. The form of the relativistic corrections adopted in the N'^LO potential of Ref. [2^ 
corresponds to the choice Ps — 1/4 and /3g = 0. The UTs driven by 5*8 and Sg also induce additional contributions to 
the current operator given by 
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(4.25) 



Evaluating matrix elements of the operators given in Eqs. (|4.22p and (|4.25p we find no contributions to the current 
density. For the charge density we obtain the following result: 
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(4.26) 



Here we have introduced fci^2 = Pi 2 +Pi,2- In addition to the constants /38,9 which parametrize the l/rriAr-dependent 
UTs and also show up in the expressions for the one-pion exchange potential, the exchange charge density in the 
above expression also depends on the arbitrary constant /33, which shows up neither in the potential nor in the 
remaining contributions to the exchange charge and current densities. We found that the corresponding UT affects 
the single-nucleon charge operator. Moreover, renormalizability of the single-nucleon charge operator enforces the 
choice /33 = 0. 
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E. Final results 



In this section we summarize the final, renormahzed expressions for the current and charge densities at order eQ, 



Ji-K = 2^ ^cx + "/tree , and piT, 
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(4.27) 



The obtained results in momentum space read 
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where the scalar functions /i(fc) are given by 
/i (fc) = 2ie— Jg , 
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and the loop function L{k) is defined in Eq. (IB. 41) . The one-pion exchange charge density has the following form: 



PlTT 



(72 ■q2 3 



92 



MS 



Si-kq2- kfj{k) + CTi • q2f8{k) 



effi 



16F^mN ql + M2 



(1 - 2^9) (t| + fi • f2) (Ti • fc<T2 • 92 



i(l + 2pg) [fi X fa] 



\ i?l • fcl0^2 ■ 92 - 0^2 • fc2(Ti • 92 ) - 2 2^, „ ^^9 0^2 ' 92 91 ' ^1 



91 + M2 

X (2^8 - l)(r| + (rl • f2)) 92 • fc + z [fi X f2f ((2^8 - 1) 9*2 • fci - (2^8 + 1) 9I • fc2 

where we have introduced 

M^ - AMI Mk) 



eg A o'l ■ 92 ^2 ■ 92 
WF^ruN {ql + M^f 

(4.30) 



/7(fc) 



h{k) 



g\ 

g\ 
'64FS 



A{k) 



fc2 



[(4M2 + fc2)^(fc) _ M^] 



(4.31) 



The loop function A{k) is defined in Eq. (|B.4 



V. SHORT-RANGE CURRENTS 



We now consider the short-range contributions. The formal structure of the currents involving the leading-order 
four-nucleon contact interactions H\^' can be decomposed into four classes 10 ... 13 as visualized in Fig. [SJ Including 
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Class 10: 




Class 11: 





Class 12: 




Class 13: 




FIG. 5: Contributions of the short-range currents. Sohd dots are the lowest-order vertices from the effective Lagrangian while 
the circle-crosses represent insertions of the electromagnetic vertices as explained in the text. For remaining notation see Fig.[T] 

the contributions induced by the UTs in Eq. (|4.2[) . we obtain the following algebraic structure: 
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^21 '/-"40 -"21 



--ff 



40 



77 — *u 77 



20 



-B^ 



-iJ, 



21 



1 7(-l).H-(l) •^^ H-(2) ^^ 7/(1) _L 1 7(-l),,H-(l) ^^ H-(2) ^^ 77(1) 



2 '-'20 '/^21 ^2 40 ^ 



i/. 



21 



P 'nH'' '-—FT 
"^20 '/-"21 77 -"' 



40 



i?2 



ff^ 



J20 '/-"21 773 -"21 '/-"40 

2-'20 '1^40 '1^21 J^3 21 



-^ /? I H-(2)„77(l) ^' 7(-l) -^^ 77(1) H-(2)„77(l) ^* 7(-l) -^^ 77(1) \ , 5 „f H-(2)„77(l) ^^ H-(l)„/(-l) 

+ P2 1 -"40 ^-"21 T;--'20 T;2"-"21 ~ -"40 ^-"21 T;J-'20 T^-"21 +P3'7 -"40 '?-"21 Ta-"21 '7-'20 



r(i)^ 7(-i)A!_ff(i) 






El' 



77(1) ^77(l)„7(-l)„ 77(2) 

-"21 773 21 '/•-'20 V-"40 



?7 + h.c. . 



Jcl3 = V 



7(0) 



7(2) 



r(i) 



•hi ;^-"4o ;^^2i ~ -^21 ;^-'^2i '7-^40 +P7?/ i -H21 ^^•^21 ^-^^40 --^^40 '?-'^2i ;p"'^2i 



r(i) -^^ 7(0)^77(2) 



(2)^77(1) -^^ 7(0) 



?7 + h.c. . (5.1) 
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We found that only the class-11 matrix elements yield non- vanishing contributions to the current density: 

dH 



Jell 



9Ai ^ 



{2^} 



■I 



21- 



k X (Tl 



T2 +ai ■ 1^2- k[Ti X T2Y 






(5.2) 



We find, however, that the resulting contribution to the current vanishes after performing antisymmetrization of the 
two-nucleon states. 

(2) 
The tree contributions emerge from gauging the subleading contact interactions in the Lagrangian Cj^-^ and the two 

new gauge-invariant terms proportional to the LECs Lx^- 



<-'C( 



e — [ri X ra] 



(C2 + 3C4 + C7) q\ - (-C2 + C4 + G-j) (cti • CT2) 91 + G-j ((72 • gi i?i + i?i • 91 CT2) 



e — — t{ [((71 + (72) X ,fi] + leLx rf ((7i - ^2) x /c + leL^ [((7i + (72) x 91] 
Id 



(5.3) 



It is reassuring to note that all the divergences of the two-pion exchange loop integrals are cancelled by the same 
redefinition of Gi that is needed to renormalize the potential [2g]. From the two LECs that are genuine to the 
current operator, only L\ gets renormalized. Employing DR and the MS-scheme, the relation between the bare and 
renormalized LEG Li has the form 



Li = Li 



2 o 4 

9a - ^9 A A 

8Fi ^^ ■ 



(5.4) 



The charge density is completely given by the class-11, 12 loop diagrams: 

dH 



pell 



„ 9a r< ^i 



f drl f ^ ^-- -- 

/ , S3 (ai-&2{k^-l^)-^i-kc72-k^ai-la2-l 



Pci2 = -e-— ^ Ct rfcTi • CT2 



dH P 



(5.5) 



In DR, the integrals entering these expressions are finite, and the result for the short-range charge density reads 



where 



Pcontact = Gt T^ (7i • A? (72 • fc /9(fc) + CTi • CT2 flo{k) 



,,,-.. ^U,,M..iMim), 



(5.6) 



fioik) 



9\ 



32FJ^ 



(M,-(4M2 + 3fc2)A(fc)) 



(5.7) 



Notice that for antisymmetric nuclear states, the two structures in Eq. (I5.6P can be combined into a single one. 
Equations (j5.3p and (j5.6p represent our final results for the short-range contributions. 



VI. COMPARISON WITH THE WORK BY PASTORE ET AL. 

We now compare our results given in the previous sections with the ones obtained by Pastore et al. [ll, 0; S] ■ Below, 
we list the (numerous) differences and, in some cases, comment on their possible origin. 



We begin with the exchange current density given in Eq. (|4.28l) . Our result for the pion loop contributions 
(i.e. terms proportional to the loop function L{k)) agrees with the one of [2| for the class-9 operator, see the 
contribution ex L{k) in f^. As Pastore et al., we also find, that the class-8 operator is cancelled by a part of the 
class-5 operator. The rest of the class-5 operator is not mentioned in 2]. As shown in section HV C[ this part 
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vanishes in our treatment due to the renormahzabihty constraint. For the seagull current, see the contribution 
(X L{k) in /3, we obtain a completely different result with even a different isospin dependence: [ti x T2]'^ as 
compared to rf 2 in [2], see Eq. (3.36) in their work. This should not come as a surprise given the fact that 
Pastore et al. did not succeed to extract the (truly) irreducible part of the amplitude for this particular topology, 
see the discussion in appendix E of i2i]. In particular, they even encountered some non-Hermitian contributions 
which then were ignored. 

• We also disagree on the tree contributions to the current density except the one ex ^21- In particular, our terms 
(X (^8,9 have a different sign. Further, we find independent contributions from both LECs ^21 and ^22, while in 
[^ they only appear in a linear combination 2d2i + ^22- Finally, Pastore et al. miss the contributions from the 
LEG dis (which accounts for the Goldberger-Treiman discrepancy) and Iq. Last but not least, we emphasize 
that our result for the current operator in Eq. (|4.28p depends on renormalized LECs dg, 9, is, 21. 22 and Iq which, 
of course, can be taken from other reactions such as e.g. pion-nucleon scattering [30[. 

• We now turn to the one-pion exchange charge density operator given in Eq. (J4.30I) . The expressions for the 
leading relativistic corrections oc l/rriN are, of course, not new and agree with the ones given in 1J|.'^ We, 
however, also obtain nonvanishing pion loop contributions to the exchange charge density, see the /r^s terms in 
Eq. (|4.30p . which are not considered in Ref. [4]. 

• Finally, our expressions for the pion loop contribution to the short-range current and charge operators also 
strongly disagree with the ones given in Refs. [3] and [J], respectively. In particular, the results obtained by 
Pastore et al. depend on both leading-order LECs Cs and Ct, while there is no dependence on Cs in our case. 
Morover, we find that the short-range pion loop contribution to the current density vanishes completely upon 
performing antisymmetrization. The origin of these discrepancies might be related to the unitary ambiguity of 
the nuclear potential and current operators. As discussed in the previous sections, we include in our derivation a 
large number of additional UTs which are possible at the given order in the chiral expansion. In particular, pion 
loop contributions to short-range current/charge operators are affected by the strong UT defined in Eq. (3.48) 
of Ref. [l2l and ^^-dependent UTs which induce additional operators listed in Eq. (|5.ip . As a consequence, 
the resulting short-range currents might be expected to be strongly scheme dependent (i.e. dependent on the a 
priori unknown angles of these additional UTs). It is the renormalizability requirement of the nuclear potentials 
and currents that provides strong constraints on the choice of the additional UTs, see the detailed discussion in 
Refs. [3, IT4I and in section |lTl and leads finally to unambiguous expressions for the (static) nuclear potentials 
and current/charge operators at the considered order. The observed differences for the short-range operators 
suggest that the results of Ref. ^] might correspond to a different choice of the additional UTs as compared 
to the one adopted in our work.'* Our findings, however, imply that such a different choice would result in 
impossibility to obtain renormalized expressions for the nuclear forces and/or the current operator. 



VII. SUMMARY 



The results of our work can be summarized as follows: 

• We applied the method of unitary transformation to work out the leading loop contributions to the one-pion ex- 
change and short-range two-nucleon electromagnetic current and charge densities. The renormalized expressions 
for the one-pion exchange charge and current operators are given for the first time. 

• We discuss in detail renormalization of the one-pion exchange contributions which provides a stringent test of our 
theoretical approach. More precisely, all emerging ultraviolet divergences have to be absorbed into redefinition 
of the low-energy constants k and di entering the Lagrangians Cl^n and >C^jv ' respectively. There is no freedom 
in this procedure as the corresponding /3-functions of all these LECs in DR are fixed and well known. We 



3 



We provide a somewhat more general result than the one of [J] by including effects due to both UTs available at this order (terms 
proportional to ;88,9)- The choice of these parameters consistent with the two-nucleon potentials of Ref. |25ll corresponds to 0g = 1/4 
and ^9 = 0. Contributions proportional to J3g are not considered in Ref. [4|. 
** For example, we could easily generate terms proportional to Cg by choosing additional UTs in a different way. 
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demonstrate, that it is indeed possible to renormalize the one-loop contributions provided one makes use of the 
freedom to employ additional unitary transformations. 

• We succeeded to obtain compact, analytical expressions for the current and charge densities both in momentum 
and coordinate spaces which can be used in future numerical calculations. 

• Finally, we provide a detailed comparison between our results and the ones obtained by Pastore et al. within a 
different framework. 

The final results of our work are summarized in Eqs. (|4.28p . (|4.30l) . (|5.3p and (15. 6p which contain the expressions 
for the one-pion exchange and short-range contributions to the two-nucleon current and charge densities at order eQ 
(leading loop order) . It would be interesting to explore effects of these novel contributions to the exchange current 
and charge densities in e.g. electron scattering off light nuclei. This work is in progress, see also t31|] for a pioneering 
calculation along this line concentrating on the two-pion exchange contributions. 
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Appendix A: Hamilton density 

In this appendix we define the expressions for the Hamilton density and currents. First let us recapitulate the 
expressions we already defined in [3| 

nt^ - ^(A^^[TX7r]7V).(7Vt[^x^]iV), 

JoV^^^ - e[7rx^]3, 
Jo2 = -e[7r X V7f]g , 

J2/°^ = e^N^cT[TXn]^N. (A.l) 

The definitions of the other parts of the Hamiltonian density and currents is given below. 



2aT- -VTrTf^f + (l-4^)f-7r(7- Vtt^ 



nl^} = —N''NN''N+—N''aN-N^aN, 
04 2 2 



N, 



[d^T: ■ a^^^2 _ (^ _ 4^) (^ . Q^^y _ (8^ _ ^)^j2 ^^ 
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0,(3) 

a/(3) 



a/(2) 

^20 



tO (1) 
•^02 



2diQ — dis 2 



9Ai 
1 






2™ 



nW^n , 



>^02 

xO (1) 
•^20 

r. (1) 

'/20 

^ (1) 
'/20 



tO (2) 
"^21 



TV 



-2ie — esafc fc • VTTa TTb , 

/ I ^ ^ ^ 7V/2 

2«ee3ab f -7:^ V-TTa k ■ VTTb + ik-^TTaWnb ] , 

epN''{deT3 + 2d7)N , 



le 

4™ TV 

ie 



■TVt 



(1 + T-"*) V + (1 + ce) (1 + r^)<T xk + 2cjaxk 



N . 



7Vt(l + r3)ViV, 



4™ AT 

-e— iV''e3Q(,T(,(d2o + ^21 rr)'^' ^^aN , 



I 
F' 



(2) 



^21 



^ (2) 
J2I 






-iidgk X VtTs - 4i(i9fc X VTTa Ta + e^abU -CT7Ta(2di6 - dis) M^ 



i{d2i ^ ) kx 



CT X V 



^22 r 



k X a 



N, 



egA 
AtunF 



N^d{Ta + 5az)T:aN . 



(A.2) 



Appendix B: Loop integrals 



The following integrals contribute to the one-pion exchange current operator at the leading one-loop order: 

2A^, 



6A,, 



(27r)3 


w; 


dH 


Z2 


(27r)3 


'-f 


d^Z 


;2 



(27r)3 Lof 
d^l ?" /'' 



(27r)3 aj+u;_(ix)+ + u;_) 



Z2'5"'' + rfc''Z3 












ri" = 6'"'B2 + k''k''B3, 



{2tt)^ uj-uj^{lu+ +w_)2 



-/'^Z^'r = {5''^k'' + 5'"'k^ + 5^''k'')Bi + k''k^k''B5 



dH 1 

(27r)3 cj>^ 

(i3; /'^/b 

(2^J4^ 



A, 



^2'5"''+yl3r/c^ 



(B.l) 



20 



where the pion energies w; and a;± are defined as 



uji = v/P + Ml , w± = '^(r±fc) +4M2 



(B.2) 



The integrals can be computed expficitly in dimensional regularization. We only give here the results for the integrals 
that are actually needed. These are: 

2fc2 \ s^L(k) 5k^ + 24M2 



3Mi 



I3 = -ICi = -2Bi = - 



127r2 727r2 

2A^ s^L{k) 5fc2 + 24M2 



3A/2 Uir^k 



21,2 



727r2A:2 



Si = -S.3 = 

^2 = - 

B5 = - 

^0 = 

A2 = 



^(/c) - 1 
87r2fc2 ' 
A^ 2i(fc) - 1 

(4M2 + fc2)L(fc) 3M2 + fc2 



87r2fc4 
A(fc) 
47r ' 
A'U + s^A{k) 



67r2fc4 



-43 = - 



Af^ - s2^(fc) 



87rfc2 



where we have introduced 



Is / s+k 



1 f k 

A[k) — — - arctan — — ;— 



^y/W+4Ml, 



(B.3) 



(B.4) 



with fc = |fc |. Further, the integral A^r is defined in Ref. [27| according to 

1 



V 1 f d'^l 
Ayr — imi — / 7 — 



-^^'-'S5^ + T<^'°<^ll> 



where the quantity L has a pole in d = 4 dimensions and is given by 



M 



d-4 r 



167r2 



1 1 



d-4 2 



(r'(l) + l + log(47r)) 



(B.5) 



(B.' 



Appendix C: Configuration-Space Expressions 



For the sake of completeness, we also give the expressions in configuration space obtained by carrying out the Fourier- 
transformation of the momentum space results 



J-(/(<fi,g2))^ / |^|^e'^"'--e^^---/(gl,g2)(2^)3<5^(gl+g2-fc 



(C.l) 



We find the following expressions: 

Ji. = {- [tI h{k) + Ti • T2 h{k)\ a2 ■ V12 



k X V12 



[fi X T2] CT2 • V12 



— fc X 



V12 X CTi 



h{k) 



k X 



V12 +ifc| X CTi 



fi{k) + ai • ( V12 + ik] ^ h{k) - la^ U{k) 



k 



( 

47rri2 



^ik-ri 



ik 



ik 



- [fi X T2] CT2 • V12 - Y I 0^1 ■ I V12 -t — 



h{k)~k{k) 



Vi; 



1/(^,^12) 



ik-R. 



(C.2) 
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where we have introduced 



f{k,ri2) 



ri2 == ri - r2 , R 



n +r2 



Vl2 = V,,, , 



(27r)3 



dx 



Jri2-{p+kx/2) 



(p2 + M2 + (1 - a;2)fcV4)2 



(C.3) 



The one-pion exchange charge density has the following form: 



PlTI 



2M^ 02 ■ Vi2r2 



(71 ■ kVi2 ■kfr{k) + V 12 -^ifsik) 



i{l + 2/39) [fi X T2]^ f CTi • fci(T2 • V12 - CT2 ■ ^2(71 ' V12 



SF^ruN 



(1 - 2/39) (r| + Ti • f2) ^1 • fc<?2 • V12 



(Ti • Vi2i?2 ■ V12 



(2^8-l)(r| + (rI-r2)) 



V12 • fc + i [fl X T2]^ I - (2^8 + 1) V12 • fc2 + (2^8 - 1 ) V12 • fci 



-M„ri2 



ik-r-i 



Stt M„ ri2 




ik-R 



(C.4) 



Finally, we also give the coordinate-space expression for the short-range currents: 



-2.g2 
ai- ka2-k -^ /g (fc) -I- cti • (72 /lo (fc) 



F2 



Pcontact — Ct T{ 
^contact - 

-|- Cj f (72 • Vi2 (7i + (7i • Vi2 (72 

+ ieL2 [ai X ^1 le*^-^(53 (^12) . 



Ak-Rx3 



S' {ri2) 



1 r^ -l3 

e — [ri X T2] 



(C.5) 



(C2 + 3C4 + C7) V12 - (-C2 + C4 + C7) (ai • a2) V12 

-e-^T^ (cti+CT2)xVi2 +ieLirf (cti - CT2) x A: 
Id L J 
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